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Abstract 

In this paper we bring into attention variable coefficient cubic-quintic nonlinear 
Schrodinger equations which admit Lie symmetry algebras of dimension four. Within 
this family, we obtain the reductions of canonical equations of nonequivalent classes 
to ordinary differential equations using tools of Lie theory. Painleve integrability of 
these reduced equations is investigated. Exact solutions through truncated Painleve 
expansions are achieved in some cases. One of these solutions, a conformal-group 
invariant one, exhibits blow-up behaviour in finite time in Lp, L^o norm and in 
distributional sense. 


1 Introduction 

In this paper we aim at studying a class of cubic-quintic nonlinear Schrodinger (CQNLS) 
equations, given in the form 

iut + Uxx + g{x, t)\u\‘^u + q{x, t)\u\'^u -|- h{x, t)u = 0 (1.1) 

in which the complex coefficients g, q and h will be assumed to have some specihc forms so 
that the equation under consideration admits a four-dimensional Lie symmetry algebra. 

The motivation for this study comes from the recent work [1] on a general class of 
cubic-quintic nonlinear Schrodinger equations given as 

iut -b f{x, t)uxx + k{x, t) Ux + g{x, t)\u\'^u + q{x, t)\u\^u + h{x, t)u = 0. (1.2) 

We performed classification of this family of equations according to Lie symmetry algebras 
they can admit. There rt is a complex-valued function, / is real-valued, k,g,q,h are 
complex-valued functions of the form k = ki{x,t) + ik 2 {x,t), g{x,t) = gi{x,t) + ig 2 {x,t), 
q{x, t) = qi{x, t) + iq 2 ix, t) and h{x, t) = hi{x, t) + ih 2 {x, t) with the assumption that < 7^0 
or g ^ 0, that is, at least one of gi,g 2 ,Qi,q 2 is different from zero. Eq. contains 

two physically important equations: cubic Schrodinger equation for k = q = 0 and quintic 
Schrodinger equation for k = g = 0 in one space dimension. 

In [T] we transformed to (jl.ip by point transformations. Therefore (II.ip appears 
as a canonical equation when classifying the family (|1.2p with respect to Lie symmetries 
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they can admit. We would like to mention main results of this paper. We showed that the 
symmetry algebra L of Eq. (|l.ll) (equivalently, of Eq. (I1.2p l is at most six-dimensional, 
that is, 1 < dimL < 6. The following results concern the canonical equation therefore 
they actually stand for the (in look) more general family (|1.2p . 

Rl. Any CQNLS equation within class (HI) having a 6-dimensional symmetry algebra is 
equivalent to the quintic constant-coefficient equation with q = qi + iq 2 , g = h = 0. 

R2. Any CQNLS equation within class dm) having a 5-dimensional symmetry algebra is 
equivalent to the cubic constant-coefficient equation with g = gi + ig 2 , q = h = 0. 

R3. The symmetry algebra of the genuine (g and q not both zero) variable coefficient 
CQNLS equation can be at most 4-dimensional. There are precisely four inequivalent 
classes of equations in this case. 

R4. None of these classes can be transformed to the standard constant-coefficient 
cubic-quintic equation with g = gi + ig 2 , q = Qi + iq 2 , h = 0. 

According to these results, when a variable coefficient CQNLS equation has a 5- or 
6-dimensional Lie symmetry algebra, it can be transformed to well-known constant 
coefficient cubic or quintic NLS equations, respectively. In the case when a variable 
coefficient CQNLS equation has a 4-dimensional symmetry algebra, we encountered four 
different classes of equations summarized in the following table. 


Table 1: Eour dimensional symmetry algebras and the coefficients in (jl.ip 


No 

Algebra 


9 

q 

h 

Li 


(51 - 


{qi + m) 

{hi -\- i/ 12 ) 

L 2 

{r,P,R,IE} 

(51 

+ m) 

{qi + iq 2 ) 

i/i2 


{P,B,D2,W} 

(51 

+ 192 ) 

{qi + 'iq 2 )t 

Ah2 
* t 

u 

{P,R,C2,W} 

{91 

+ m) 

{qi+iq2){l+t‘^) 

Q(i+t^) 


The basis elements of subalgebras are given as follows: 

T = dt, p = d^, W = d^, B = td^ + ^xd^, (1.3) 

Di = xdx + 2tdt - ^Pdp, D 2 = ]^xda, -h tdt - (1-4) 

Ci= xtdx + t^dt-^tpdp + ^x^di_o-, C2 = xtdx + {I+ t^)dt - tpdp + ^x’^dui- (1.5) 

Here u : ^ C is expressed in terms of the modulus and the phase 

u{x,t)= p{x,t)P^^^'^\ (1.6) 

Equations o with coefficients from Table 1 are representatives of class of equations 
which admit non-isomorphic four-dimensional symmetry algebras. If we are to say that the 
symmetry algebras do not extend to five- or six-dimensional ones at all, we have to put the 
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conditions {g,h) 7^ (0,0) for Li, {g,h2) 7^ (0,0) or (g',/12) / (0,0) for L2, {g,h2) 7^ (0, |) or 
{Q,h2) 7^ (0,{0, i}) for L3 and (g,/i2) / (0,0) for L4 [I]. Let us note that Li is non-solvable, 
L2 is nilpotent, L3 and L4 are solvable algebras. Besides, Li is decomposable whereas 
the others are not. Eq. with coefficients from Table 1 cannot be transformed to a 

constant-coefficient equation, therefore they will be the main subject of this article. 

Before moving to obtaining reduced equations for each representative equation for 
algebras Li, L2, L^, L4, let us mention another noteworthy point for which the results of 
this paper may be of use. Suppose one is interested in finding radial solutions to 

iut + Au + g{x,t)\u\‘^u + q{x,t)\u\^u + h{x,t)u = 0 , (x, t) € x M. ( 1 . 7 ) 

Radial solutions of (ll. 7 p satisfy 

iut + Uxx - Ux + g{x,t)\u\'^u + q{x,t)\u\^u + h{x,t)u = 0, x>0 (1.8) 

X 

where x plays the role of the radial coordinate, u = x^^^u transforms (|l.ll) to (|1.8p (up to 
a change in the coefficients). Thus using the results of this paper one can easily produce 
the corresponding results for the family (jl.Sp and hence for the multidimensional one (jl. 7 p . 

2 Reductions Through One-Dimensional Subalgebras 

Now we proceed to find the reduced ODEs for any representative equation (|l.ll) for 
Li, L2, L3, L4^ with coefficients from Table 1 . Complex equation (II. 1 |) has a four-dimensional 
solution space with two dependent and two independent variables. In order to reduce the 
solution space to a three-dimensional one with two dependent and one independent variable 
(i.e., a system of ODEs), we only need one-dimensional subalgebras of Li, L2, L3, L4^. Thus 
we restrict ourselves to subalgebras of each symmetry algebra having generic orbits of 
codimension 3 in the solution space. 

The classification of one-dimensional subalgebras under the action of the group of inner 
automorphisms of the four-dimensional symmetry groups is a standard one. We do not 
provide the calculations leading to the conjugacy inequivalent list of subalgebras. The 
classification method can be found for example in laEnn. We can refer to [ 5 ] for the 
classification of subalgebras of real three- and four-dimensional Lie algebras, from which 
the content of the following Table 2 is a direct consequence. 


Table 2 : One-dim. subalgebras of 4 -dim. algebras under the adjoint action of the full 
symmetry group 


Algebra 


Subalgebra 

a, 6, c G M, e = =Fl 

Li 

Li.i = {T + Ci+ aW} 

L1.2 = {Di + bW} 

L1.3 = {T + cIE} 

L2 

L2.1 = {P} 

L2.4 = {IE} 

L2.2 = {T + aW} 

L2.3 = {B + bT} 

L3 

L3.1 = {Pj 

L3.4 = {D2 + aW} 

L3.2 = {B} 

L3.5 = {IE} 

L3.3 = {P + eP} 

u 

II 

L4.2 = {C2 + alE} 

II 

CO 
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In view of the formulation (HSI) we express dm) as a system of two real equations for 
the modulus p and the phase oj as follows 

- pujt +Pxx - +giP^ + qiP^+ hip = 0, (2.1a) 

Pt + “^Px^x + P^^XX + g2P^ + Q2P^ + h2P = 0. (2.1b) 

2.1 Non-solvable algebra Li = {T, Hi, Ci, IT} 

Let us note that the algebra has the direct sum structure Li = sl(2,M) © {W} and the 
representative equation of the algebra is 

1 1 

iut + Uxx + {gi + ig 2 )-\u\'^u + {qi + iq 2 )\u\'^u + {hi + ih 2 )^u = 0 (2.2) 

X 

with the real constants gi,g 2 ,qi,q 2 ,hi,h 2 . 


Subalgebra Li,i Invariance under the subalgebra Li.i implies that the solution has the 
form 


/ ^ M{C) 

u{x,t) = —exp 
^X 


i ( a arctan t + 


x'^t 


4(1+ f2) 


+ P 




X 


1 + P 


(2.3) 


and the reduced system of equations satisfied by M(^) and P{C) are 


/ 3 _ c2 . 

AfM" - 4fMP'^ + -aC + hijM + giM^ + qiM^ = 0, (2.4a) 

4f{M'^P'y + h2M'^ + 52 + 52 M® = 0. (2.4b) 


We first need to decouple these equations to solve for the functions M and P. It is seen 
that an integral of ()2.4bp can be obtained for two different cases of the constants. 

(i) 52 = g '2 = 0 , /l 2 7 ^ 0 

It can be shown that the system (|2.4I) amounts to integrating a third order nonlinear 
ordinary differential equation from (I2.4ap : 


1 , 


y'y'” - + zy'y” + ^ 

2 ^ 2^2 

where the functions M, P of 


1 / 3-^2 


-aC+hi)y'2-^y'3 


h2 -^2 

are related to y(^) by the relations 

Am = 2ifiP', = 


(2.5) 

( 2 . 6 ) 


(ii) 52 = 52 = ^2 = 0 

In this case we can easily decouple the reduced system of equations. Integration of 
()2.4bp gives 

M^P' = C, P{0 = I C = const. (2.7) 

and from (j2.4aD we obtain the equation for M 


4^2 4^2 


1 a 

+ *16 + 4^" 


3 + 4hi 
16^2 


M + C^M 


-3 


( 2 . 8 ) 
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For the remaining subalgebras, we only state the results obtained through similar operations. 


Subalgebra Li 2 


u[x, t) = —— exp 


i{b\nx + P{0) 


, ^ = 


t 


(2.9) 


results in the system 

4^2M" - - 46)MP' + (f + hi - b‘^)M + giM^ + qiM^ = 0, (2.10a) 

^2 4P') ' + + g2M^ + q2M^ = 0. (2.10b) 


(i) 52 = 92 = 0, h2 / 0 

1 , 


Y'y'" — —Y"^ + —Y'Y” + 
2 ^ 

-^iY + Cf = 0, 


3 + 4hi _ J_W/2 _ _9]_y,3 , Qi^y>4 

8^2 4^ 32/ 2/12 2/i2 


1 


= P =- - - Inf- 

h2 00 “. 


e ^ . h2 fY + C 


^Y 


rdC, 


(2.11a) 

(2.11b) 


(ii) g2 = q2 = h2 = 0 


M" = C^M-^ + {—- ^ ^ (2.12a) 

^8^ 16^2 64^ 4,^2 4^2 ’ V y 


e b 


dC 


P= h _ _in^-C / — 

82^ / M 2 


(2.12b) 


Subalgebra Li 3 

results in the system 


t(x,/) = M(x) exp i(ct + P(x)) 


M" - MP'2 + (^ - c)M + —+ giM® = 0, 


x\M^P'^ + /i2M^ + 52 xM^ + g2X^M'’ = 0. 


(i) 52 = 92 = 0 , /i 2 7 ^ 0 

y'y'" - -y"2 + -Y'Y" + 2(^ - c)y'2 - ^y'3 + ^^y'^ 

2 3/ V 3/ ' ^2 ^2 

2 / 7 ^ 

-^(y + c)2 = o, 

= "^2 J da:. 

(ii) g2=q2 = h2=0 

M" = C^M-^ +{c-^']m - —M^ - qiM^, 

\ J X 

f C 

= Tridx- 


(2.13) 

(2.14a) 

(2.14b) 


(2.15a) 

(2.15b) 
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(2.16a) 

(2.16b) 



2.2 Nilpotent algebra L 2 = {T, P, B,W} 

The canonical equation has the form 

iut + Uxx + (ffi + ig-2)\u?‘u + {qi + iq2)\u\'^u + ih2U = 0. (2-17) 


Subalgebra L2.1 The invariant solution u{x,t) = M(t)exp leads to the first-order 

system 

P' = M' + h2M + g2M^ + = 0 (2.18) 

which is straightforward to integrate upon given values of the constants. 


Subalgebra L 2,2 u{x,t) = M{x) exp i{at + T’(x)) leads to the system 


M" - MP'^ -aM + giM-^ + qiM^ = 0, 

{M'^P'y + /laM^ -h g2M* + gsM® = 0. 


(i) 52 = 92 = 0 , /l 2 / 0 

Y'Y'" - -y"2 - 2ay'2 - ^y'3 + ^y'4 _ 2hl{Y + Cf = 0, 
2 h2 h2 

M{x) = {-y/') ' ’ = -^2 I dx. 

(ii) g2 = q2 = h2=0 

M" = + aM - giM^ - qiM^ 

= / ^dx. 


(2.19a) 

(2.19b) 


( 2 . 20 a) 

( 2 . 20 b) 


( 2 . 21 a) 

( 2 . 21 b) 


Subalgebra L2.3 

(A) The case 6 7 ^ 0. The invariant solution 


u{x,t) = M(^) exp 
gives the system 


1 1 3 r. 

—xt - + P 

2b 662 


e 



M" - MP'2 - -4jM + = 0, 

2^4 

b^iM'^P')' + -h 52 -h 52M® = 0. 


(i) 52 = g '2 = 0 , 6-2 / 0 


r/\r/// ^ ^ X ^/2 251 


2qibP 


64' 


h2 


hi 


262 

y’y'" - ^y"^ - + ^4^y'^ - -^(y + c)2 = o, 


62 . 


M(e) = {—-Y'fl\ 
6-2 




Y + C 
Y' 


dp 


( 2 . 22 ) 

(2.23a) 

(2.23b) 

(2.24a) 

(2.24b) 
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(ii) g2 = q2 = h2=Q 


M" = C^M-^ + -IjM - - |iM^ 


P{x) = I 


264 


62 - 


62 - 


C 

M2 




(B) The case 6 = 0. 


2 

u(x,t) = M(t)exp i(^^ + P{t)J 


(2.25a) 

(2.25b) 

(2.26) 


is the form of the group-invariant solution and the reduced system of equations is 

P’ = 51 m2 + q^M\ M' + {h2 + ^)M + 52M3 ^ ^ q (2.27) 

This system is readily solved by standard methods. 

2.3 Solvable algebra L 3 = {P, B, D 2 , W} 

We note that the algebra is not decomposable and the representative equation from Table 
1 is 

iut + Uxx + (51 + W 2 )\u\^u + (51 -I- iq 2 )t\u\^u + i—u = 0. (2.28) 


Subalgebra L3.1 u{x,t) = M(t)exp (zP(t)) leads to the system 

P' = 5iM2 ^ M' + ^M + g2M^ + q2tM^ = 0. 


Subalgebra L3.2 u{x,t) = M(t)exp 


4t 


gives the reduction 


(2.29) 


P' = g^M^ + q^tM^, M' + 


26-2 + 1 


2t 


M + g2M^ + q2tM^ = 0. 


Subalgebra L3.3 u{x,t) = M(t)exp 


i( I 


m) 


gives 


(2.30) 


P' = g^M^ + 5itM4, M'+ {— + 


M + g2M^ + q2tM^ = 0. (2.31) 


t 2{l + et)_ 

Subalgebra L3.4 A group-invariant solution invariant under L 3,4 will be of the form 


1 r 1 

tt(x,t) = —M(^) exp z( 2 alnx-|-P(^)) , ^ = — 


(2.32) 


where M, P will be solutions of the system 


4^2 M" _ 4fMP'‘^ + ^(^ - 8a)MP' - 2^M' -h 2(1 - 2a^)M -h 51 M^ -h yM® = 0, (2.33a) 
4.f{M‘^Py - 2^M^P' - - 8a)MM' -h ( 62 C - 6a)M‘^ -h 52 ^^ -F ^M® = 0. (2.33b) 
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()2.33bp is satisfied if 52 = <?2 = 0, /12 = 1/4, P{^) = - — aln^ + Pq and (I2.33ap reduces to 

8 

1 . a 1 1 


M" = —M' + (—-- -)M - —M^ - —M^ 

2e ^4^ 2^2 64^ 4^2' 

2.4 Solvable algebra L 4 = {P, B, C 2 , W} 

The last canonical equation under investigation will be 


2 \| |4 , ■ t P h2 


iUt + Uxx + (51 + m) \u\‘^u + {qi + iq2){l + t^)\u\'^u + i 

Subalgebra L 4 ,i u{x,t) = M{t)exp + P{t)) 


2 ( 1 +t 2 ) 

gives the system 


u = 0. 


(2.34) 


(2.35) 


r2 I „ n I nw , l + /l2l + 2t^ 


P' = g^PP + qi{l+P)M^, M'+ 


Subalgebra L 4 2 


M+ g2M''^+ q2{l + t^)M^ = 0. (2.36) 


2 t(l + t 2 ) 


u{x, t) = 


1 


ViTP 

results in the system 


M (^)exp 


i ( a arctan t + 


x'^t 


4(1+ t2) 


+ P 




X 


VTTP 


M" - MP'^ - (a + j)M + giM^ + qiM^ = 0, 

(M^P')' + y + g2M^ + q2M^ = 0 . 


(2.37) 

(2.38a) 

(2.38b) 


(i) 52 = (72 = 0 , /i 2 / 0 


1 ^2 


o ^^W2 4^1 .3 ^ 8^1 .4 hi 


Y'Y'” - -y "2 - (— + 2 a)y'^ - ^y'^ + -'p^{y ^cf = o, (2.39a) 

2 2 /i 2 ^ 

/ 2 ^ 

M (0 = { y 




h2 fY + C 
Y' 


di. 


(ii) g2 = q2 = h2=0 


M" = - 51- qiM^ 


C 




(2.39b) 

(2.40a) 

(2.40b) 
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3 Analysis of the Reduced Equations 


In this section we will analyze the reduced equations obtained in the previous section. 
We wonder whether any of those reduced equations are integrable, so as to obtain exact 
solutions for them, therefore to the PDEs they are originating from. The main tool for 
the search of integrability is going to be the Painleve property of differential equations and 
related techniques based on this concept. We would like to mention a few lines about the 
methods we will utilize first. 

A differential equation is said to have the Painleve property if all its solutions do not 
include movable singularities other than poles. Painleve himself classified all second order 
ordinary differential equations in the complex plane which have this property laE] and his 
work was completed by Gambler [8] . There are fifty canonical second order ODEs having this 
property. All but six of them are integrable in terms of known functions. The six canonical 
differential equations, named Painleve DEs I-VI, define new transcendental functions. A 
complete list of these equations can be found in [9]. This is the case for second-order 
ODEs, but not the end of the story. The classification was continued for higher order ODEs 
afterwards (see for example [iniE]). When we have a second- or higher-order nonlinear 
equation in hand integrate, or to look solutions for, checking any correspondence with these 
Painleve-type equations is the very first station to stop by. This is going to be our first 
approach in analyzing the reduced equations. 

In order to check whether a given partial differential equation has the Painleve property, 
a necessary condition was proposed in [12]. According to the so-called WTC test, the 
unknown function u{x, t) of a PDE (system) is expanded into a power series in the form of 
a Laurent series 

OO 

u{x,t) = ttj(x, t) d>(x, (3.1) 

j=0 

Here ^{x,t) = 0 is an arbitrary singularity manifold (a curve in this case). The 
power a is called the leading order. If a assumes (negative) integer values, the Laurent 
series form excludes singularities other than poles. In the case of a PDE system, upon 
plugging these form of expansions for dependent variables one obtains a system of the form 
Q(j)U(j) = F(j). It is not possible to determine coefficients Uj for which detQ(j) = 0. 
They remain free, if this system is compatible for those indices jres, called as resonances. 
The partial differential equation is said to pass the Painleve test if there arises the correct 
number of arbitrary functions as required by the Cauchy-Kovalewski theorem given by the 
expansion coefficients, and ^{x,t) should be one of the arbitrary functions. 

If a differential equation passes the P-test, it has a good chance of integrability and one 
can try transformation to the representative equations of the canonical classes mentioned 
above, besides making use of the results available in literature for them. If the consistency 
conditions at the resonance indices are not satisfied, the PDE does not pass the Painleve 
test, therefore it does not have the Painleve property. In this case, if one truncates the 
P-series at a finite term as 


N 

u{x,t) =Uj{x,t) ^{x,t)‘^~^^ (3.2) 

j=0 
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and try to satisfy the consistency conditions for special cases of functions Uj and (which 
failed to remain free), (13.211 may turn out to be an exact analytical solution for the PDE. By 
this way one can also obtain Backlund transformations relating solutions of the equation. 
We refer the interested reader to la o 15] and m, among many other works on 
Schrodinger and other type equations. 

After this overview of the methods we shall make use of, we can turn back to have a 
look at the material we obtained in the previous section. Table [3] is a guide to the equations 
obtained as reductions for the canonical PDEs with symmetry algebras Li, L 2 , L^, L 4 . In 
what follows, we list our findings about those reduced equations. 


Table 3: A summary of the reduced equations 


Order 

Equation Number 

1 

(I2.18D, (I2.27D, (I2.29D, (I2.3UD, (I2.31D, (12.361) 

2 

(12.81), (l2.12aD. (l2.16aD, (I2.21al), (l2.25aD, (I2.34D, (I2.40al) 

3 

(12.51). (I2.11al), (l2.15aD. (I2.2()al), (I2.24ap. (l2.39aD 


• The first order decoupled systems can be solved by elementary methods, so we do not 
present any solution for them here. 

• Among the second-order equations, it was seen that (I2.21al) passes the Painleve test 
after the transformation M = y/U without any restriction on the constants. Let us 
note that, (I2.21ap is a reduction for the PDE 

iut + Uxx + gi\u\‘^u + = ^- (3.3) 

Ref. [I7| is devoted to the search of exact solutions of (|3.3p . and ()2.21ap is obtained 
as a direct reduction, after all being completely analyzed to obtain elliptic solutions. 
We refer the interested reader to this work. Note that since (I2.21ap does not include 
the independent variable explicitly, it can be integrated to a first-order equation by 
the transformation M' = W{M). 

• The second order equation (I2.16al) passes the Painleve test after the transformation 
M = y/U if gi = hi = 0 and the previous argument is also valid for this equation. 

• Second-order reductions other than (|2.16ap and (|2.21ap fail to pass the Painleve test, 
therefore they give pretty less expectation for integrability. 

• Prom the third-order equations, (|2.20ap does not include the independent variable; this 
means it can be reduced to a second-order equation by a transformation Y' = W{Y). 

• All the third-order equations pass the Painleve test if the variable Y is replaced as 
Y —)• 1/F{^) or y —?• 1/F'{^). After this change of variable, the third order equations 
convert to an equation in F which has a leading behaviour with a = —\ or a = —2 
with several branches successfully satisfying the conditions of the test. 

As we have mentioned, among seven second order equations we obtained, only (|2.16ap and 
(I2.21al) passed the Painleve test under some restrictions on the parameters they include. 
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However, under these conditions, the PDE (13.311 they originated from turned out to be 
one which was already analyzed in literature. Therefore, our first approach for analyzing 
the reduced equations, which is to compare them with the representatives of P-integrable 
equations of second order did not give any new result to us. At this moment we take another 
approach. We look for solutions through truncated P-expansions for our equations. Section 
3.1 is an effort to do so for second-order reduced equations. Section 3.2 tries the same 
method not for the reduced equations but by directly attacking the canonical PDEs for the 
algebras given in Table 1. 

3.1 Exact solutions to second-order reduced equations throngh truncated 
expansions 

All of the second-order equations listed in Table [3] has a leading order a = —1/2. This 
suggests that the transformation M = \/U will transform those equations to equations with 
dependent variable U, for which the P-series expansions will have a leading-order a = — 1. 
Eor the moment let us denote the independent variables for second order equations in Table 
[3] by X. The Painleve series expansion for solutions U{x) will have the form 

OO 

U{x) = Uj{x)^-^+^{x). (3.4) 

j=0 

We truncate this P-series expansion at the first term; that is, we propose a solution to our 
second-order ODEs in the form 

C(.) = (3.5) 

When we plug this expression in the DE for U, terms which can be collected according to 
the powers 4>-^, j = —4, —3, —2, —1,0 appear in the equation. By putting the condition that 
the coefficients of those <h-^s vanish, we try to determine the form of the functions Uq{x) and 
<h(x) and yet determine U{x). For the seven second order ODEs in the Tabled we have 
been successful in two cases. 


Subalgebra Li, 3 . We start with Eq. (|2.16ap . including the constants gi,qi, hi of the original 
PDE, the constant c arising from the group action and the integration constant C obtained 
during the decoupling. Putting M = i/U in (I2.16ap gives 

„ 1/9 / hi ^ 2qi 9 n 2C^ 

U" = —U’^ + 2 c - 4 H _ 2 q^u^ + 3,5 

2U x^ X U 


We impose (|3.5I) in (13.6p and require that coefficients of ^-^s arising vanish. First of all, by 
the coefficient of it is necessary that qi is a negative number, say qi = —/?^. Further 
calculations put several constraints on the constants included in the equation. We must 
have c = C = 0. We see that it is possible to solve the overdetermined system of equations 
for Uq and ‘h in two different cases. If 


gi 


27/3 

V3’ 


qi = -/ 3 ^ 


hi 


2 _ 2! 

2 4 ’ 


7 7^ -1 


11 









we obtain 


<i>(x) = + k2 


Uo{x) = 


V^ki{l + 7) 
2/3 


(3.7) 


where [3 and 7 are fixed by the given values of gi, qi and ki,k 2 are arbitrary constants. 
Similarly, if 

-2/3 2 , 3 

51 = qi = -p , hi = -- 


we find 


[/o(x) = 


^/3ki 1 


2/3 X 

We state the final solution as follows. 


<l>(x) = kilnx + k2- 


Proposition 1 The canonical equation 

I 5 i+/ 52 | |2 I / . N| ,4 , hi+ih2 

iUt+Uxx3 -|w| U + [qi + iq 2 )\u\ u3 - ^ 


u = 0 


for the subalgebra Li = sl(2,M) (BW is solved by the following functions, 
(i)- If 9i = qi = -/3^ = 7 - X' 7 7^ -1; 52 = 92 = = 0 ; 


u{x,t) = Cix) = A 


-11/2 


kix^+'f + k2 

— -2/3 „ _ a2 h. — 3 


pPo 


2 _ V^kiil + 7) 


= 


2/3 


(ii). Ifgi = qi = -/3^ hi = c?2 = 92 = /12 = 0 : 

Ae^^° 


u{x,t) = 


2 y/ski 


[x{ki Inx + ^2)]^/^ ’ 


= 


2/3 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


Remark 1. In fact, (|2.16al) is of Painleve type if gi = hi = 0. The above solution does 
not have such a requirement, therefore is a solution in non-Painleve case. 


Subalgebra L2.2- Next we give the analysis for (|2.21al) . It includes the constants gi and 
qi originating from the reduced PDE, whereas a is a free reduction parameter. Setting 
M = y/U results in 

1 9/^2 
u" =—U'^+ 2aU - 2giU^ - 2qiU^ + —. (3.12) 

Again it is necessary that qi = —/3^ < 0. We find 

o2 7^ 

qi = -P , ® Y 


Uo{x) = ^>(x) = + /c2. (3.13) 


We proceed by plugging (|3.5p in (13.120 . 
that for 

27/3 

the functions Uq, will be 


^/3fc1 O' 
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Proposition 2 The canonical equation 

i'^t “ 1 “ ^xx + (91 + i92)\u\^u + {qi + iq2)\u\‘^u + ih2U = 0 
for the subalgebra L 2 is solved by 


u{x, t) = A 




kie^^ + k2 


1/2 




\/3fci7 

~w~ 


where gi = qi = g 2 = q 2 = h 2 = 0. 


(3.14) 


(3.15) 


Remark 2. One needs to keep in mind that (|2.21ap is a Painleve-type equation without 
any restriction on the constants. 

We tried the same procedure for all other second-order equations, arriving at failure. 


3.2 An exact solution by truncation in Painleve series 

The method we are going to use is truncation in Painleve series again. Unlike the previous 
subsection, this time we are focusing on a canonical PDE itself, instead of its reductions. 
We would like to note an exact solution of the variable coefficient canonical equation 

, , fl'l+* 52 | |2 , / , . M |4 hi+ih2 to ^c:\ 

lUt + Uxx^ - \u\ u + iqi + iq 2 )\u\ u-|- 5 -m = 0 (3.16) 

X x^ 

of the symmetry algebra Li. For the case qi = q 2 = 0, inspired by m, we obtained 
interesting exact solutions in |15j . giving rise to blow-up analysis studied in [El- 

In fact, when we make a leading order analysis to (j3.16p . it is immediate to see that 
it cannot have the Painleve property as the balancing of the leading orders is not possible 
because of the quintic term. Even the leading order a cannot be determined. However, we 
keep on tracking the route we followed in the absence of the quintic term technically (see 
m for details), we write (I3.16P with its complex conjugate v as the system 

, 9 i + i 92 2 , t , • \ 3 2 , to -1^ \ 

tut + Uxx^ - u V + [qi + zq 2 )u V H- 5 - u = 0, (3.17a) 

X x^ 

, 9i ~ i92 2 . t • \ 3 2 I ~ *^2 „ /„ 1.71 N 

-ivt + Vxx^ - V u + [qi - iq 2 )v u H-- v = u (3.17b) 

X x^ 

and propose a solution of the form 

u{x,t) = uo{x,t) {x,t), v{x,t) = vo{x,t)^~^~^^^{x,t). (3.18) 


The system (I3.17P has no leading order. Those powers are proposed in regard 

to cubic case studied in [15]. When we substitute (I3.18P in (I3.17p . terms with coefficients 

system. We shall note here first two of them: 


{qi + iq 2 )ulvl^ ^ + 

{qi - iq2)ulvl<^-^^^^ + 


91 + m 

X 

' 91 - m 

X 


'Uq'Uq “t“ (2 H- 

■uqVq + (2 — 3i(5 




4> 


— 3 — 25 


$ 


—3-l-i(5 


(3.19a) 

(3.19b) 
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In a usual truncation procedure one proceeds with requiring that coefficients of and 

vanish separately. However, it is not possible to have that coefficients of 
vanish. We collected these terms together in the form 


{qi + iq 2 )^ ^ UqVq + ttoUo + (2 + 3i6 - 

{qi - iq 2 )^~^ ulvl + ——— uqVo + (2 - 3i6 



uod>-3-*‘5 = 0, 





(3.20a) 

(3.20b) 


and required that their sum vanish. These equations imply that we must have 


(|)2 

UqVq = A —, d'(x, t) = (t)o{t) X 

X 


B 


(3.21) 


where 

^ 3.5ffi + ((5^ - 2)g2 „2 ^ +gi^ 

3 (I(?i + (<52 - 2 ) q 2 ’ <52-2 


(3.22) 


If we proceed and consider that coefficients of <I)“2 =f*< 5 vanish, we see that we 

need to have B = 2/3, hi = 5/36, /i 2 = 0 and we are able to determine the form of uq, vq 
and successfully as follows. 


<h(x,t) 

uo{x,t) 

vo{x,t) 


^-1/2 ^2/3 


(cit + C2)2/3 ’ 

3 , 1/6 

{Cit + C2)2/3 """P 
3 , 1/6 

{cit + C2)2/3 """P 


ClX 


_ 4(cit + C 2 )J ’ 

Cix 2 
— I 


4(cit + C 2 ) 


(3.23a) 

(3.23b) 

(3.23c) 


We need to re-formulate (I3.22h using the condition B = 2/3. We skip the details here and 
sum up the main result. 


Proposition 3 The function 


u{x, t) = 



ClX 


4(cit + C 2 ) 


-din 


/cx2/3 


{cit + C2)2/3 


(3.24) 


solves (13.1611 for hi = 5/36, /i 2 = 0 where 5 


— ^(<72 + 922 I) and A is a root of the equation 


32 - 9A{-4gi + A{-4qi + {g 2 + q 2 Af)) = 0. (3.25) 


Remark 3. Proposition 1 also presents a solution to Eq. (|3.16p . (|3.10p seems to contain 
(j3.24p when the magnitude is considered. However, Proposition 3 does not have the 
condition 32 = g 2 = 0, therefore (I3.24p is valid for a more general choice of the parameters. 
Furthermore, the solution which is the main theme of [TB], in which the cubic equation is 
considered, coincides with (13.241) in special cases. 

We tried the same approach for the canonical equations of L 2 > E 3 , L 4 but this ended up 
without any success. 
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4 Blow-up in Lp, L^o Spaces and in Distributional Sense 


The machinery of this section will follow that of [16] . The main consideration here is going 
to be the solution (I3.1U|) of Proposition 1 found for the canonical equation (13.9p . Let us 
stress what kind of importance (I3.9p has by the following proposition. 

Proposition 4 Any PDE of the form dni) (or, equivalently which admits 

a four-dimensional Lie symmetry algebra containing sl(2,M) as a subalgebra can be 
transformed by point transformations to the canonical form (1331). The symmetry algebra of 
(1331) is spanned by the vector fields T, Di,Ci,W noted in ([L3ll-(fT3D. 

These vector helds produce point transformations of translation in time, scaling, conformal 
and gauge transformations, respectively. By exponentiating the generators T, Di and Ci 
and combining them we find the SL(2,M) group action on the solutions to (|3.9I) . We state 
this result in the following Proposition. 


Proposition 5 Given a solution Us{x,t) of (13.91) . 


u{x, t) 


, ,1/9 i . ( X c + dt\ 

+ bt) e Us I -—,- 1 


(4.1) 


is also a solution to ()3.9I) for ad — be = 1. 

Now we take the solution Us{x,t) to (13.9p as C{x) defined by ()3.9p : Us{x,t) = C{x) of (13.9p . 
If we use Proposition 5, we have that 

u{x, t) = (a + bt)~^^‘^ e*4(a+6t) ^ f —^ (4.2) 

\a -\-bt J 


also solves (j3.9l) . Suppose o > 0, 6 < 0 and define e = a-\-bt = b{t — T), where T = — ^ > 0. 
Let us note that t —>■ T~ if and only if e —)• O'*'. This enables us to simplify (14.2p in the form 

C,e{x) = C (^) • (4.3) 

Let us write down C,e{x) explicitly: 


C,(x) = A 






(4.4) 


Some facts for later convenience: This expression has no discontinuity for x > 0 if j, ki, k 2 
and e are positive. But x'^ or may not be dehned for x < 0. However, this can be 
prevented by putting the condition that 7 = 2 ^ — 1 >0 where ni, n 2 are natural numbers 
relatively prime. In this case, the solution for x > 0 and the solution for x < 0 are the same 
in modulus but they have a phase difference of With this condition, becomes an 
even function. On the other hand, if we do not want to put this condition on 7 , then we 
may consider the function 


Ce(x) 


Ce(x), X > 0 
0, X < 0. 


(4.5) 
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Let us call Ce(a;) or Ce(®) as u^{x). In case any question arises on whether is indeed 
a solution in regard to differentiability issues at x = 0: By a direct calculation we see 
that u^{x) is twice differentiable with respect to x at x = 0 for e > 0 if 7 > 4 (the value 
of the second derivative being zero). Therefore, Ue{x) is a classical solution to ()3.9p on 
—00 < X < 00 , t < T for 7 > 4. Observe that 


/ 


CO poo 

Uf;{x)\dx = cr \ue{x)\dx 

Jo 


(4.6) 


where d = 2 if ttg = Ce and d = 1 if = Ce- 


4.1 Blow-up in Lpy 
Theorem 1 For any T > 0 there exists a solution u{x, t) to Eq. (|3.9p such that 


lim ||u(x,t)|L = 00 , p>2 

t^T- 

where ||tt(x,t)||p ~ J W\^dx'^ 


t^T- 
i/p 


(4.7) 


Proof. Let T > 0 be given. We will take u{x,t) as Ue(x) we have defined. We are free to 
choose numbers a, b such that a > 0, b < 0 and T = — — . We can say that 


lim ||u(x,t)|L = lim ||ue(x)|L. 
t^T- e^0+ 


Let us calculate Lp(M) norm of Ue{x). 

\w{xWp=aw r 

Jo 

If we substitute x = ey, we get 


poo 

x'^ 

Jo 

iklX^+'T + k2€^+'y\ 


lP/2 


dx. 


(4.8) 


(4.9) 


I i Mip r 

I 


e 2 Jo 
a\A\P ^ 

P-2 

e 2 


y 


.7 


lP/2 


+ k2- 

7 


y 


+ k2- 


dy 

p/2 roo , 

dy + 


lP/2 


I lkiy^+'y + k2- 


dy 


(4.10) 


Clearly, the first integral is finite due to the continuity of the integrand and the second 
integral is finite iff p > 2. Therefore for p > 2 we get 


r II t ill 
hm \\u,{x)\\p = — 


>0+ 


g 2p 


( r 

\ 1 

[Jo 

-kiy'^y + k2- 


lP/2 

dy ] = 00 . (4.11) 


This completes the proof of Theorem 1. 
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4.2 Blow-up in Loo(K) 

The following theorem concerns Ue{x) and shows that it blows up in Loo(]R) norm. 

Theorem 2 Given any T > 0, there exists a solution to (IH.Qp such that 

lim ||w(x,t)||oo = oo. (4-12) 

t^T- 

Here ||n(x,t)||oo = esssup^gjo,^) \u{x,t)\, t < T. 

Proof. We again concentrate on the solution Ue{x), for which we similarly write 


lim ||u(x,t)||oo = lim ||we(x)||oo- (4-13) 

t^T- e^0+ 


By the definition of ne(x) we can restrict our attention to the interval [0,oo). Since |ue(x)| 
is continuous on [0,oo), |u£(0)| =0 and lima;^oo \ueix)\ = 0, there exists xq G (0,oo) such 
that 


|^^e(3;)l|oo = ess sup \Uf:{x)\ = max |ue(x)| = |j4| 
3;S[0,oo) x£[ 0 ,oo) 


feixj+^ + k2e^+'r 


1/2 


(4.14) 


We find that this maximum occurs at 

= G15) 

and 

7 

Kn r> ^ 1 + ")' 

||u,(x)||oo = Kl= ^ _^ -. (4.16) 

A:/+^A:2^+^(1 + 7) 

Therefore we arrive at 

lim ||ue(x)||oo = oo. (4.17) 

6^.0+ 

Before the conformal transformation, the solution (j3.10p is a stationary solution. The 
conformal transformation gives the time-dependent solution (14.4p . The stationary solution 
is obtained by putting e = 1. By the action of the conformal transformation the maximum 

of the stationary solution at the point x = becomes a moving one as it is in 

(I4.15P and the maximum point approaches the origin as e ^ O"*" t —)• T~, whilst the 

maximum value inreases unuboundedly. This is the case if ttg = Ce- If = Ce there are two 
peaks at approaching origin and blowing up. 


4.3 Blow-up in the space of generalized functions 

We denote the space of infinitely differentiable functions with compact support by D and 
D' is the dual of D, the space of generalized functions. 
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Theorem 3 When p > 2 


P-2 


e 2 \ue{x)\^—> K5{x) as e ^ 0"^ 


(4.18) 


in D' where K = 
the origin. 


poo 

/ ^ 1 ^ 1 ^ 

40 


kiy^+'y + k2 


p/2 


dy and 5{x) denotes the Dirac distribution at 


Proof. We need to show that 


P-2 , 


lim / €2 \u^{x)\^(f{x)dx = Kip{0), V(/9 E D. 

<5^0+ 


(4.19) 


Let us denote the integral on the left by /(e). 


/ OO _2 _2 /*C 

\ue{x)\^ip{x)dx = e^/ 
-CO J 0 


+ /c 2 e^+">' 


p/2 


ip{x)dx. (4.20) 


Plugging in X = ey gives 


poo 

/(e) = / 

4o 


kiy^+'y + k2 


p/2 


ip{€y)dy. 


(4.21) 


Let (t\AY’ 


p/2 


kiy^^'^ + k2 

|/,(y)| < C,a\A\P 


ip{ey) = feiy). Since cp is of compact support, we have 


y 


kiy^^'^ + k2 


p/2 


lim /,(y) = a\A\P 

e^0+ 


y 


kiy^^"< + k2 

By Lebesgue’s dominated convergence theorem we find that 

poo poo 

lim / fe{y)dy = /^(O) / a\AY 

4o Jo 


£Li([0,oo)) if p > 2, 

p/2 

yp(O). 


kiy^^'^ + k2 


p/2 


dy = //<y9(0), 


which implies 


P-2 


e 2 \ue{x)\P —)• K5{x) 


(4.22) 

(4.23) 

(4.24) 

(4.25) 


as e —O’*". 

Let us finally stress that a similar machinery for blow-up in distributional sense is carried 
out in HSUS] for Davey-Stewartson type equations. 


5 Conclusion 

In this paper we analyzed a class of variable coefficient cubic-quintic nonlinear Schrodinger 
equations. We concentrated on representative equations of families with Lie symmetry 
algebras of dimension four. What is special about these equations is that they cannot 
be converted to standard constant-coefficient cubic-quintic NLS equation. They do have 
variable coefficients. In a methodological way, we obtained several different reductions of 
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these PDEs to ODEs. Among them for the ones being Painleve-integrable, we referred to 
results existing in literature. We have seen that many of the reduced ODEs inherited the 
non-integrable character of the relevant PDE, leaving us with less chance for presenting exact 
analytical solutions. However, we still have been successful for doing so, using truncated 
expansions of Painleve series. 

The truncated expansions were tried in obtaining exact solutions both for each of the 
four canonical PDEs themselves and for their seven second order reductions. The attack 
was fruitful for two of the second-order reductions by producing three different solutions, 
and for one of the PDEs itself by producing one exact solution. 

That being not the end of the story, we used one of these exact solutions to expose the 
blow-up character of solutions of one of the canonical PDEs, which includes inhomogeneous 
nonlinearities and an inhomogeneous potential. The blow-up character is due to conformal 
invariance of the equation, and occurs in Lp norm of the solution with p > 2, in Lqo norm 
and in distributional sense. All solutions we have considered are to be taken in the classical 
sense. 

We would like to finalize by noting how the results of this paper may be of (further) 
use. We have presented some exact solutions for variable-coefficient PDEs. They will be 
useful when the PDE is to model a real phenomenon. If not at all, these solutions can be 
used for the convergence test of numerical algorithms of PDEs. We have obtained many 
reduced ODEs of second and third order which did not have the Painleve property therefore 
not giving us the chance to integrate them by using the known results in literature for 
P-integrable equations. Being born on earth from a well-known equation of mathematical 
physics, these equations may be analyzed for possible classical and non-classical symmetries 
(e.g. A-symmetries). What is more; these reduced ODEs may also be useful in cases where 
the numerical schemes for the PDE they were obtained require high computational cost or 
produce divergences, as there are many sophisticated numerical algorithms for the treatment 
of an ODE [20]. Finally, to our best knowledge, the search for formation of blow-up in 
suitable Sobolev spaces remains as an open problem. 
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